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Abstract: In this paper, the effects of periodic partial harvesting of a continuously grown 
crop on augmentative biological control are analyzed. Partial harvesting can remove a 
proportion of both pests and biological control agents, so its influence on the control efficiency 
cannot be a priori neglected. An impulsive model consisting of a general predator-prey model 
in ODE, augmented by a discrete component to depict releases of biological control agents 
and the periodic partial harvesting is used. The periods are taken as integer multiples of 
each other. A stability condition for pest eradication is expressed as the minimal value of 
the budget per unit time to spend on predators. We consider the partial harvesting period 
to be fixed by both the plant's physiology and market forces so that the only manipulated 
variable is the release period. It is shown that varying the release period with respect to 
the harvest period influences the minimal budget value when the former is carried out more 
often than the latter and has no effect when releases take place as often as or less frequently 
than the partial harvests. 

Key-words: Predator-prey dynamics; partial harvest; inundative biological control; im- 
pulsive control; stability; minimal budget 
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Effets de la recolte partielle des cultures sur la lutte 

biologique 

Resume : Ce document montre les effets de la recolte partielle periodique sur la lutte 
biologique inondative dans le cas de la protection d'une plantation a croissance continue. 
Comme la recolte partielle est susceptible de prelever une partie des ravageurs ainsi que des 
auxiliaires de lutte biologique, on ne peut ignorer son influence sur son efncacite. 

Le modele etudie consiste en un modele proie-predateur en EDO classique decrivant 
l'interaction biologique entre les deux populations, auquel s'ajoute une partie discrete re- 
prsentant les phenomenes, intrinsquement discrets, lies a la recolte et a la lutte biologique. 
Les parametres de recolte partielle sont supposes fixes par la physiologie de la plante ainsi 
que la demande economique; la periode de lacher des auxiliaires et le budget investi sont 
done les seuls parametres modifiables du systeme. L'exigence de stabilite de l'etat du sys- 
teme ou les ravageurs sont absents nous donne un budget minimal d'auxiliaires a utiliser par 
unite de temps qui peut etre fonction de la periode des lachers. Nous demontrons que la 
periode des lachers influe sur le budget minimal quand elle a lieu plus souvent que la recolte, 
mais n'a aucun effet quand elle a lieu moins frequemment. 

Mots-cles : Dynamiques proie-predateur; recolte partielle; lutte biologique inondative; 
commande impulsive; stabilite; budget minimal 
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1 Introduction 

Biological control is the reduction of pest populations to harmless levels through the release 
of their natural enemies. The latter can include both parasitic and predatory species, which 
are deployed at selected locations throughout the crop and, wherever possible, to specific 
parts of individual plants where the pest is likely to attack. Successful control projects in 
the field have involved the use of only one predatory species such as in [3[ [9] , as well as 
more complex biodiverse schemes such as those suggested by [l4t [SJ [2J [17] and the references 
therein. The target pest species and the setting, i.e. where the crop is grown, usually 
determines the type of control required, namely whether pest eradication is necessary or 
not. For an exhaustive list of definitions and applications, we refer the reader to [5} 116). 

In this report, we consider the protection of continuously grown crops which have zero 
tolerance to pest invasions. There are two aspects in this type of culture. 

• Firstly, inundative control which is a prophylactic method of pest control yields the 
most satisfatory results when implemented (see [3 [20l [6] for theoretical/simulatory 
studies and [H [HI El H] for real life experiments) . A calculated number of predators 
are repeatedly injected into the ecosystem, independently of the detection of pest 
insects in the greenhouse. Such populations are not allowed to thrive and consist only 
of individuals whose main source of subsistence is the pest insect, in the absence of 
which, they (the predator insects) rapidly die out. The frequency of the releases and 
the number of predators injected each time ensures that a minimal 'sentry' population 
is present to reduce the damage caused by the pests on their attack. 

• Secondly, over their growing period, these crops are partially harvested on a regular 
basis. Since it is known that harvests are likely to influence, even counterintuitively, 
predator-prey dynamics |18[ I15j. it has to be taken into account in the formulation of 
the problem. 

We consider the simplest ditrophic case whereby one predatory species is used to eradi- 
cate a pest population. Our model consists of ODEs augmented by a discrete component to 
incorporate the effect of partial harvest and releases that by their very nature are discrete 
phenomena. This is a classical formulation that is used widely in the literature where impul- 
sive dynamics are studied. Examples are |12l I15j in the context of agricultural ecosystems, 
|19j in epidemiology, [10] in pulsed chemotherapy to cite some. Few papers in the literature 
on impulsive crop protection however seem to focus on stability of the pest-free state: yet 
this is of practical importance especially for high valued crop cultures. 

In our work, we attempt to give an economic dimension to the solution of our problem 
by defining the releases in terms of the number of predators to invest in over a budget 
period. Using Floquet Theory as presented by [19], we are able to express the stability 
condition as the minimal number of predators per budget period required to drive the pests 
to zero at a given release frequency. [13] showed how this number varied with the release 
period chosen. The worst case scenario of pest attack occurring at an intermediate stage 
between two predator releases was considered and the optimal release policy which would 
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guarantee the most efficient protection against surges in the pest population was calculated. 
In particular, it is shown that the higher the frequencies of predator release, the smaller the 
time interval over which the pest population was above a threshold commonly referred to 
as the Economic Injury Level [21] - and hence the lower the damage incurred by the crops. 

In line with the work of [13J, we investigate how the frequency of releases is to be varied 
with respect to the (fixed) harvesting frequency to minimise the minimal budget value. 
We consider the harvest period as a reference since it is set by market constraints. The 
effect of partial harvesting is similar to that of pesticide usage proposed by [12] in their 
Integrated Pest Management (IPM) strategy. Our model departs from the latter's in three 
ways. Firstly, both the predator and pest populations are subjected to partial harvesting 
when this occurs. Secondly, hypotheses made on the functions governing the population 
changes are weak and can encompass most of the density-dependent functions proposed in 
the literature. Finally, one period is taken as the integer multiple of other. This feature 
is key in solving for the stability condition to obtain the minimal budget value. The case 
where the frequencies are not the same is included. 

It is shown that for a given harvest period, when releases take place less often or as often 
as harvests, the minimal budget is at a calculated value which is independent of release 
period. However, when releases take place more often than harvests, the minimal budget 
required always exceeds this value. This result runs counter with that obtained by [T3] : 
merging the two seems to indicate that the harvest frequency is a threshold that should not 
be exceeded when releasing predators for efficient biological control. 

In the first section of this article, the system model is presented. The mathematical 
analysis of the system's stability and the formulation of the stability condition in terms 
of the minimal budget are presented in the next section. A brief interpretation of the 
mathematical results follows. Finally, we conclude with a discussion on their implications. 

2 Model description 

The model we present consists of a continuous part to depict the predator-prey interaction. 
We consider the case at the onset of pest invasion where the crop - the pest food supply - 
is in abundance. Because of this, at this stage, it is sufficient to model only the pest x and 
predator y species. 

f(x) - g(x)y 
h(x)y - dy 

(1 - a x )x{nT h ) VneN (1) 

(1 — a y )y(nTh) + 5 (nTh mod T r ) uT r Vn 6 N 
(1 — S (mT r mod Th) a y )y(mT r ) + uT r Vra G N 

The first two equations govern the intrinsic predator-prey interaction occurring in the 
system. The three ones depict the impulsive phenomena that we consider with harvest 
taking place at nTh and releases at mT r . 



y 

< x(nT+) = 

y(nT+) = 

. y(mT+) = 
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In the continuous part, the functions discussed are not specified so they are representative 
of as many systems as possible. Only the following hypotheses are made. 



Hypothesis 1 Let f(x), g(x) and h(x) be locally Lipschitz continuous in R + such that 

• /(0) = 

• 0(0) = 0, g'(0) > and g(x) > Vx > 

• h(0) = and h(x) > Va; > 



f(x) is the growth velocity or feeding input of the pests. It represents the growth function 
of the pest species and in our model, it also encompasses any non-predatory losses of the pest 
population (e.g. logistic growth). We assume that the predator population is never large 
enough for intra-predator interaction to take place so the functional and numerical responses 
can be expressed solely in terms of the prey numbers, i.e. as g(x) and h(x) respectively. 

We assume that pest growth rate, the functional and numerical responses are all nil when 
the ecosystem is pest-free. 

The functional response is increasing for small pest population levels. We also consider 
that, in the presence of pests, predation always takes place with a negative impact on x 
(g(x) > 0) and a positive impact on y (h(x) > 0). Note that conditions can be induced as 
much by the predator insect foraging abilities per se as they can be facilitated by placing 
the predator insects at known locations on the plant where the pests are most likely to 
attack. In classical density dependent models, g(x) is bounded or linear, so that is 
always bounded. The boundedness of ^1 means that there is no value of x where the 

■> g(x) 

pest growth f(x) overwhelmingly dominates the predation g(x), which would render the 
biological control impossible. 

Partial crop harvests and predator releases occur respectively every Th and T r . a x and 
a y represent the respective proportions of the prey and predator populations affected at 
each harvest. These parameters are allowed be different since in reality, it is very likely that 
each species tends to occupy different parts of the plant. We also assume that the insects are 
uniformly distributed throughout our plantation so that the effect of partial harvesting is 
directly correlated with the number of plants harvested. We assume linear maturation of the 
crop so the proportion of crops harvested each time and hence insects removed is considered 
as fixed. The <5-function is defined thus to identify instants of simultaneous partial harvest 
and predator release. 



Finally, we presume that we have a fixed budget of predators over a designated time 
period that is distributed evenly among the releases that are carried out. \i refers to the 
total number of predators purchased per time unit. Expressing T r in the same units as the 
budget period gives the control fiT r as the number of predators released every T r . 





(2) 
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3 Mathematical analysis 

In our analysis, we restrict ourselves to the case where either one of the periods (release 
or partial harvests) is the integer multiple of the other. Note however that the model {T]) 
formalism is more general. We study the system in the absence of pests, i.e. when x = 0. In 
addition of being invariant, it is the target state of our system. The stability of the system 
around that state is therefore of interest. Our analysis takes place separately for the case 
when releases are more frequent than harvests, and when they are less frequent. 

We show that in the absence of pests at the initial time, the predator population converges 
towards a positive periodic solution. We then demonstrate that when preys are present at 
the initial time, convergence of the predator population also takes place to that same periodic 
solution, while the preys go extinct provided some condition on the parameters is verified. 



3.1 Pest-free stability analysis 
Releases more frequent than harvests 

Proposition 1 Let = kT r where k 6 N* and Hypotheses [7] be satisfied. Then, in the 
absence of pests, model (Qp possesses a globally stable periodic solution 



(x vh (t),y ph (t)) 



l m ;, " J-i 



0)y * e -<i(i mod T„) + ^-dit mod TV) ^ 



(3) 



3=0 



whe 



1 - (1 - a v )e- dT » 



fJ.Tr 



(4) 



Proof: When Th = kT r , in the absence of pests and using Hypotheses [U the system is 
simplified to 



x =0 
y = -dy 

x(mT+) = (1 — S (to mod k) a x )x{mT r ) 
2/(toT+) = (1 — 8 (to mod k) a y )y(mT r ) + fj,T r 
Vm e N 



(5) 



The pest population stays nil since in the absence of pests, their population does not 
change either. The solution 

x p h{t) = 



is trivial. 
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On the other hand, the predator population will vary according to the number of preda- 
tors manually injected into the system and, since the population is non-zero, according to 
the partial harvest effect. The absence of their source of food will cause an exponential decay 
of the population. We demonstrate that these forces will provoke the predator population 
to reach a periodic pattern of period equal to Th, which we shall refer to as the period of 
reference. The instant following a coinciding partial harvest and release is taken as the point 
of reference. 

To prove Proposition [TJ we first show by induction that the predator population right 
after a release can be expressed in terms of the point of reference as follows 

i— l 

y(nT h + iT+) = y(nT+)e- idT * + fiT r £ (6) 

3=0 

where i G [0,1,..., (k- 1)] 

It is seen that |(6]) is valid for i — since it is equal to 

-l 

y(nT+) = y(nT+)e° + fiT^e'^ = y(nT+) 

3=0 

Now suppose that holds for i = q where q £ [0,1, ... ,k — 2], i.e. 

9-1 

y(nT h +qT+) = y(nT+)e"^ + +^T r e~ jdT - 

3=0 

We will now show that is valid for % = q + 1. We calculate y(nTh + (q + 1)T+) from 
y{nTh + qT+) using ©, then substituting (|3.ip into ([7]) as follows 

y(nT h + {q + l)T+) = y{nT h + qT+) e - dT ^ + fiT r 

y{nT+)e-* dT r + ^T r ^ ] e~ dT ^ + fiT r 

3=0 ) 

= y(nT+)e-^ +1 ^ + pT r £ e"^ + /xT r 
= y(nT+)e-^ dT " + fiT r £ 

3=0 

so that |(6]) holds true for i € [0, 1, . . . , k — 1]. 

To evaluate the evolution of y according to the period of reference Th , we need to calculate 
the value of y((n + 1)T^), which is equivalent to y{nTh + kT+), in terms of y(nT^) . At 
this point however, partial harvesting takes place before predator release; so we first express 
it in terms of y(nTh + (k — 1)T+) then expand the expression using |(6|) as follows 



(7) 
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y((n + l)T+) = y(nT h + (k-l)T+)e- dT "{l-a y )+fiT, 



fe-2 



,{nT+)e- d ^-^ + \iT r £ e - dT "(l - a y ) + M T r 

3=0 J 

fc-1 

y{nT+)e- dT * + fxT r (l - a y ) £ e^' dT " + M T r 



j'=i 
fc-i 



2/(nT+)e- d ^+/xT r ( (1-%)^ 

j=0 



Note that the summation term can also be evaluated so the sequence is expressible as 

/ 1 - e~ dTh \ 

y((n + l)T+) = y(nT+)e- dT -+^T r \^l-a y ) Y ——^r+a y j (8) 

In this linear dynamical system, the coefficient of y(nT£), e~ dTh is less than one in magni- 
tude, so the sequence will converge to a limit, the equilibrium of ([8]). This equilibrium yields 
(J4j) and the convergence of y(t) to a periodic solution y p h(t) based on y*. 

Now that we have established the existence of the periodic solution y p h{t), we seek to 
formulate it. We focus on a reference period over nT/, < t < (n+l)Th during which y p h(t) is 
piecewise continuous, with the continuous components separated by predator releases. The 
continuous intervals are defined over nTh+iT r < t < nTh + {i + l)T r where i € [0, 1, ... , k— 1]. 
For a given value of t, the value of i is easily identified as being i = [ 1 m ^ d Th J . The value 
of y P h(t) is then of the form 

y P h(t) = y ph (nT h + iT+)e- d{ - 1 mod T "> 

and, from ^ with y{nT^) = y* , we have that 

i-i 

y ph (nT h + iT+) - y*e^ + fj,T r £ e^ dT " 

3=0 



so that 



V P h{t) = [y*e- idT r+(iT r J2 



jdT r I -d(t mod T r ) 



e 



3=0 



y * e ~d(t mod T h ) + ^Y r e~ d{ ~ 1 mod Tr) ^2 eT 2dTr 

3=0 



y * e ~d(t mod T h ) + ^J^g-dC* mod Tr) e _idTr 
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This is of the same form as proposed in (J3j) , thereby completing our proof. □ 
The form of the y p h function is illustrated on Figure [TJ 




T r 2T r T h = 3T r 4T r 5T r 2T h 7T r 8T r 3T h 

t 



Figure 1: Form of the periodic solution y p h{t) in the case where k — 3. Releases of predators 
are apparent at every mT r instant, while the cumulative effect of harvest and release leads 
to an apparent smaller release at every nT/j instant. Between those instants, the population 
decays exponentially since it has no prey to feed on. 



Releases less frequent than harvests 

When harvesting is more frequent than the release of predators, we have a similar result 
about the existence of a periodic solution. 

Proposition 2 Let T r — kTh where k € N* and Hypotheses [7] be satisfied. Then, in the 
absence of pests, model (Qp possesses a globally stable periodic solution 

(x pr (t),y pr (t))=(0,y*e- d ( t mod T ">(1 - % ) L ^T^ J ) (9) 
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where 

fiT r 

y = i n \k -dT ( 10 ) 

1 — (1 — oi y ) K e a±r 

Proof: When T r = kT^, in the absence of pests and using Hypotheses [TJ the system is 
simplified to 

x =0 
y = -dy 

x{nT+) = (l-a x )x{nT h ) (11) 
y(nT^) = (1 - a y )y(nT h ) + fj,T r 8 (n mod k) 
Vu eN 

As previously explained, x pr {t) is solved for trivially as being 

x pf — 

We prove that the predator population again reaches a periodic solution. This time, 
however, the period of reference is T r . The point of reference is the instant after a coinciding 
harvest and release. We show by induction that the population after a harvest can be 
expressed as 

y(mT r + %T+) = y(mT+)e^ dTh (1 - a y Y (12) 

where i G [0,1,..., (k- 1)]. 

It is seen that lfT2l is valid for i — since it resumes to 



y(mT+) = y{mT+)e°{l~a y ) 
Suppose (fl2|) holds for i — q where q G [0, 1, ... k — 2], i.e. 

y(mT r + qT+) = zj(mT+) e -'? dTfc (l - a y )" (13) 



We will now show that (fl2|) is valid for i = q+1. We calculate the value of y when i = q + 1 
in terms of y(nTh + qT+), knowing from y — — dy in fTTj) it will be an exponential decay 
with the added component for the harvest. We then substitute lfT3|) in and obtain 



y(mT r + (q + 1)T+) = y(mT r + gT+^-^Cl - Oy) 

= (y{mT+)e-i dTh {l-a y ) q )e- dTh {l-a y ) (14) 

= y(nT+)e- rf («+ 1 ) r "(l " 

This is clearly the same form given from the expression in lfT2|) . thereby validating it. 

y* is given as the fixed point of the sequence representing post-release instants. Therefore, 
using (fT2l) for i — k and model ifTTj) . we next calculate y((m + 1)T+) as 

y((m + l)T+) = y(mT r + kT+) 

= y(mT+)e- kdTh (1 - a a ) fc + (15) 
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In this linear dynamical system, the coefficient of y(mT+), e~ dTr (l — a y ) k is less than 
one in magnitude, which confirms the existence of the fixed point y* to which the sequence 
converges. This equilibrium yields (fT0|) and the convergence of y(t) to a periodic solution 

y pr (t). 

Now that we have established the existence of the periodic solution y pr (t), we seek to 
formulate it. We focus on a reference period over mT r < t < (m + l)T r during which y pr (t) is 
piecewise continuous, with the continuous components separated by harvests. The intervals 
of continuity span mT r + iTh < t < mT r + (i + l)Th where i € [0, 1, . . . , k — 1]. For a given 
value of t, the value of i is easily identified as being i — [ 1 Tr j . The value of y pr (t) is 
then of the form 

y pr (t) = y pr (mT r + iT+)e- d ^ mod T ^ 
and, from ifT^I with j/(mT+) = y* , we have that 

y pr (mT r + iT+) = y*e~ tdTh (1 - a v ) 1 

so that 

VpAt) = (y*e- ldTh {l-a y ) l )e- d ^ modT ^ 
= y*e~ d(t modTr )(l-a/ "" J 

which is exactly the expression given in ([9]) and completes the proof. □ 



The form of the y pr function is illustrated on Figure [H 



3.2 Global stability analysis 

Since we will study the convergence of the solutions to (0, y p (t)) (where the p subscript 
stands as well for ph or pr), it will be convenient to describe the system in terms of the 
deviation coordinates with respect to the reference periodic solution: 

x(t) = x(t) — x p (t) 
y(t) = y(t)-y p (t) 



This yields 



and 



x = f(x)~g(x)y 

= f(x)-g(x)(y + y p (t)) (16) 



h(x)y -dy - h(x p )y p + dy p 

h(x)(y + y p (t)) - dy (17) 
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3T h 4T^ 5T h 



Figure 2: Form of the periodic solution y pr (t) in the case where k = 3. Harvests are 
apparent at every nTj, instant, while the release of predators dominates the harvest at every 
mT r instant. Between those instants, the population decays exponentially since it has no 
prey to feed on. 



The impulsive effects on x are obviously unchanged compared to those on x. On the 
other hand, the release effects on y disappear in y; indeed, we have 

y(mT+) = y(mT+) - y p (mT+) - y(mT r ) + ^T r - (y p (mT r ) + M T r ) = y(mT r ) 

The harvesting impulses are preserved in the expression of y 

y( nT h ) = y( nT h ) - V P ( nT h ) = ( 1 ~ a y )y(nT h ) - (1 - a y )y p (nT h ) = (1 - a y )y(nT h ) 

In the sequel, we will perform a global and a local stability analysis. For the latter, we 
will need the computation of the linear approximation of the deviation system around the 
periodic solution (0,y p (t)): 

i - (f(0)-g'(0)y p (t))x 

i = h'(0)y p (t)x-dy 1 } 
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Releases more frequent than harvests 

We will first prove our result in the case where releases take place more often than har- 
vests. We obtain two different constraints for the Local Asymptotic Stability (LAS) and 
Global Asymptotic Stability (GAS) of the periodic solution in system ([I]). The latter is 
obviously stronger than the former, but is sufficient in the case where pests outbreaks do 
not immediately take large proportions. 

In order to state the following theorem, we first need to define the function 

t± h (S,r) = d[S+ — — 



a y (l-e- dT h) 
l-(l-a y )e- dT h 



e~ fc 
k[ 1-e k 



This function is increasing in S and r because the sign of the partial derivatives is determined 
by the sign of the last factor, which can be shown to be positive. Indeed, this factor is positive 
when 



a y (1 - e- dT ") 



_ilh 
e ~* 



l-{l-a y )e-^) U(!_ e -^ 
and we have o^e" - ^ < a y and 1 — (1 — a y ) e~ dTh > a y , so that 



< 1 



a v e k 

< 1 



1 - (1 - a y )e 



-dT h 



Also, k(l- e-T^J > (1 - e' dTh ) since both sides of the inequality have the same value in 
Th = 0, and 

A M 1 --*))- A < 1 - , "*> 



which shows that fj, (S,r) is increasing in S and r. 

Theorem 1 When Th = kT r with k £ N*, the solution (x(t),y(t)) = (0,y p h(t)) of (Qp is 
LAS iff 

M> ^(Si' 5 ' (0) ) (19) 



and is GAS if 

yx>Q g[X) x > X 
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Proof: We start with the proof of global convergence under condition ((20| . In this proof, 
we will first show that x goes to zero, from which we will derive that y goes to also (so 
that y(t) converges to y p h(t)). 

Let the initial condition for system i|16p - fT7j) be (io,2/o) at time to = + , that is after 
the harvest and the predator release that take place at the initial time. Analyzing (flT)) and 
noting that y p h{t) + y — y(t) > 0, we have 

y > -dy 

so that y(t) > min(0, y )e~ dt . 

In order to analyze the x equation, we define the function 

G(x) = [ X -±-ds (21) 
Jx 9{s) 

which can easily be seen to be an increasing function of x since g(s) > 0. Since we also have 
that 

g(s) < (sup x>0 ^r-) s, it is straightforward that lim G(x) = — oo. In order to show the ex- 

x— >0 

tinction of the pests we will then prove that G{x) goes to — oo as t goes to infinity. Therefore, 
we write the G dynamics: 

dG{x) = 1 ~ 
dt g{^) 

= w\ ~ y - yph( f ) 

< ^-min(0,y )e- dt -y ph (t) 

We will now consider the evolution of G between two successive harvests, that is the evolution 
of G between the times nT^ and (n + l)Th for a given n: 



r{n+l)T h 

G(x((n + l)T h )) < G(x(nT+j) + 



-min(0,y )e s - y P h(s) 



ds 



Since no impulse is present inside the integral, we can drop the + superscript in its lower 
extremity. 

We will now analyze how the harvest that takes place at time (n + 1)T^ impacts G. We 
have 



G(*((n + 1)T+)) = C n+1)Tit) jk ds 



_ rx((n+l)T h ) i , f 2((n+X)T+) i , 
- Jxo g(s) aS ■+■ Jx((n+l)T h ) g(s) aS 

< G(x(nT+)) + f^ 1)T * {ggjJ - min(0,y o )e- ds - y ph (s) 



ds 



+ Jx{{n+l)T h ) g(s) aS 

(22) 
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The last term represents the influence of harvest on G and can easily be approximated be- 
cause 

SU Pz>0 ^) 



x((n + l)T h ) > x((n + 1)T+) = (1 - a x )x((n + l)T h ). Denoting S g = sup„ >n *M and 



r g = sup^^o we have 

l-(l-a x )x((n+l)T h ) l 
Jx((n+l)T h ) 



9{s) 



ds < 



(l-a x )i((n+l)T h ) 



x((n+l)T h ) 



1 h^l-a*) 
— as = 

r n s r„ 



(23) 



Introducing ([23)1 into (|22|) then yields a bound on the application between successive 
moments after harvest. 



G(x((n+1)T+)) < G(x(nT+))- 



(n+l)T h 



- min(0,j/ )e ds - y P h{s) 



nT h 



ds+ 



ln(l - a x ) 

r 9 
(24) 

We can now evaluate an upper-bound for G at any time t > 0. Defining / as the integer 



WW) 



part of j^, we have: 



ds 



, ln(l-« x ) 



G(x(t))-G(x ) < /o[58S -min(0,y o )e- ds -y p/l ( S )_ 

< /„' [S s - min(0, y Q )e- ds - y ph (s)] ds + iMl^l 

= - $1 min(0, y )e- ds ds + f* Th [S g - y P h{s)] ds + I J^ h [S g - y P h{s)] ds 

I ; ln(l-a x ) 

= ^§^l(e' dt - 1) + j* Th [S B - y ph (s)] ds + I ^ [S g - y ph (s)] ds 

I ? ln(l-^) 



The first two terms are bounded (the first one is obvious and the second one is upper-bounded 
by S g Th). We then have to analyze the third one, which has been obtained through the 
periodicity of y p h(t) and the fourth in order to know if G(x(t)) goes to — oo when t goes to 
infinity. In fact, it suffices to have 

^ [S g -y ph (s)]ds+ Hl ' ax) <0 

r g 

to achieve this. It is more cleanly rewritten in the form 



Th mm ^ c t , m (! ~ ax) 
y p h(t)dt > S g T h H 



(25) 



In order to obtain (|2"U)l . we are now left with the computation of f^ h y p h(t)dt, which is 
detailed in Proposition [3] of the Appendix: 



y P h{t)dt 



^T h 



1 - 



a y (1 - e 



-dT h 



_ilh 



1 - (1 - OL y )e- dT * 



k ( 1 - e"T 



(26) 



RR n° 6284 



16 



S. Nundloll, L. Mailleret & F. Grognard 



Introducing lf26|) into ([25| then yields ([20| . which shows that this last condition is suffi- 
cient for having x going to as t goes to oo. 

Since x goes to zero, there exists a finite time tf after which h{x) < | for all times. 
Therefore, after this time, we have 

V = h{x){y P h{t) +y) - dy < h(x)y ph (t) - ^y 

We have seen that h(x)y p h{t) goes to zero as t goes to infinity; so does also y. 

In order to have the global asymptotic stability, we are only left with the local asymptotic 
stability to prove. In order to do that, we only have to consider the discrete system that 
maps the state at time nT^ onto the state at time (n + l)Tj[ with respect to the linear 
equation ([TH)) and the discrete part. After some computations, we obtain: 

*)((n + l)T+)= B (|)(«T+) (27) 

where 



B 



i (i 



C»+l)T h 



Note that % is a term that we do not use in our analysis, therefore is not expressed. 
Indeed, since the matrix is triangular, it is stable if \B\i\ < 1, i.e. 

(n+1)T \, dr , /'(Q)r* + Mi-«») m , 

^ y ph dr > m (28) 

Similarly to what was done earlier, it can be shown that ([281 is equivalent to tfT~Q|) . so 
that the necessary and sufficient condition for local stability is proven. 

It is directly seen that f28|) is satisfied when 11251) is because fj, (S,r) is increasing in S 
and r and we have 

/'(0) ,. f(x) . f(x) g(x) g(x) 

— — = hm — — — < sup — — and g (0) = lim < sup (29) 

9 (0) x ^ g(x) x >q g{x) x ^ x x > x 

This completes the proof of global stability, since we have shown global convergence and 
local stability when lj20j) is satisfied. 

□ 



Releases less frequent than harvests 

If we now consider the case where predators releases take place less often than harvests, we 
also obtain global and local stability results based on the following function 

Br (S,r)=diS+ j 1 _ e _ dTh 
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which is increasing in S and r since the last fraction is positive and a x < 1. 

Theorem 2 When T r = kTh with k 6 N*, the solution (x(t),y(t)) — (0,y pr (t)) of (Qp is 
LAS iff 

M>M,f£S,$W) (30) 



-'As'(o) 

and is GAS if 



i f( x ) 9( x )\ 
M>M, sup— - T , sup : (31) 



r \x>o g{x) ' x>0 X 

Proof: This proof does not depart very much from the one of Theorem [TJ The only- 
difference is that the reference period is now T r . We use the same function G[x) as in (|2"Tj) 
and an analysis identical to the one of the previous theorem leads to 



G(x(mT r + {l + 1)T+)) < G{x{mT r + IT+)) 



ds + 



ln(l— qj 



which is exactly (|24| since it is depicting the behaviour of the model between two harvesting 
instants. 

Extending this to the whole T r interval, we obtain 



G(x((m + 1)T+)) < GiximT+^ + J^^^^-miniO^e-^-y^s) 



ds 



We now see that this expression is identical to (|24| with the exception of the presence of a 
k factor and the expression of y pr (t), which comes from instead of ([3]). 
Condition lj25j) then becomes 

f y pr {t)dt > S g T r + fc M 1 -^) ^ 
JO r g 

and the computation of J^ r y pr (t)dt with y pr (t) as in ^ yields: 

Tr uT 1 - e~ dTh 

^ ( ^ = T hi^ ( 33 ) 

This leads to conditon ([3lj) (the computation of lf33|) is detailed in Proposition [4] in the ap- 
pendix). Global convergence of (x, y) to (0, 0) is then concluded by using the same argument 
as in the proof of Theorem Q] to show the convergence of y to 0. 

The local stability condition lj30j) then directly arises from the analysis of the stability 
of the discrete linearized system that maps y(mT+) onto y{{m + 1)T+). 

*)((m + l)T+)=B(j)(mT+) (34) 
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where 

__/(!- a x fe^ )Tr fVi-SWtor*- \ 

~ \ ef^t 1)Tr (1 - a y fe- d &l l)Tr ^ ) ^ 

Again, since the system matrix is lower triangular, for stability we simply require that 
\B\x\ < 1 and I-B22I < 1- The latter yields a trivial condition, so we calculate the former: 

(m+1)T " . , f(0)T r + k\n(l-a x ) 

VprdT > — (36) 

which directly leads to condition f30|) for local stability and the proof of global stability is 
also complete because /i is increasing in S and r and f29|) is still satisfied. □ 



Comment 

As we have seen, when the condition ([20|) or IpTTj) is satisfied, the extinction of the pests is 
GAS. When the local condition lfl9l) or f30|) is not verified, the extinction of the pests is not 
stable and a bifurcation analysis similar to what is done in [lOl [T5] would show the presence 
of a limit cycle when \x is close to the limit. When fi satisfies condition (fT9|) or ([301 only, the 
pests extinction is locally stable and we cannot rule out that it is globally stable (since our 
global condition is only sufficient). Such a budget has the advantage of being smaller than 
the one that guarantees global stability. It allows for good control of limited pest invasions; 
however the culture is at risk of being destroyed by a large pest outbreak. 

Since, in both cases, the conditions for local and global stability are identical up to 
two different parameters, any analysis of the consequences of one of those conditions will 
immediately translate to the other. The interpretation of conditions l|19p - l(20|) and f30|) - ([3Tj) 
will be given in the next section. 



4 Interpretation of results 

It is easy to see that \i is independent of T r . The influence of T r on \i is trickier to identify 
so we shall analyse it mathematically first. We then present graphically the variation of 
both [i^ and [i h with respect to T r for a typical set of parameter values, and attenpt to give 
a practical interpretation of these results. 

4.1 Mathematical analysis 

We first need to note that when S+ '"^"^ < 0, for any of the local or global condition, the 
condition is trivially verified. Indeed, it implies simply that no biological control is needed 
for exterminating the pests; in fact, the partial harvesting is effective enough for this purpose 
(as a x is large enough). We now evaluate how the release frequency influences the minimal 
budget when this condition is not trivial. 
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We have already seen that fi is independent of T r . We will now study the latter's 
influence on fi . 

Theorem 3 Let T h = kT r where keW. 

The minimal budget is monotonically decreasing with respect to the release period T r for 
non-negative values of fih, i-e. 

w r < < 37 > 



rp OfA 

Proof: Knowing that T r is equal to it is possible to identify the sign of -ff- noting 
that 

d t h _ % dk _ % 



dT r dk dT r dk \ T h 

So 



dT r J ° V 9k 

fi, is expressed as the product of two distinctive parts, one of which is independent of k and 
which, for the non-trivial stability condition, is positive, 

s+ Hl - a * )d >o 

rT h 

where S and r are the parameters required for the local and global conditions, as defined 
previously. 

The second part is viewed as a composite function of k so that f38|) can be evaluated as 



where a{k) = ( k ^_l-*T h/k \ ) ■ Then, we get 



tLh\ _ l-(l~a y )e- dT h da 



fda s 



Sgn {dk) (40) 
/ e- dT "/ fc fdT\ _ _ dTh/k 

^ SgI1 U 2 (l-e^ Th/fc ) 2 V k 
-sgn(fce- dT "/ fc + dT h - k) 



RR n° 6284 



20 



S. Nundloll, L. Mailleret & F. Grognard 



Since 

d 



( ke ~dT h /k + dn _ k] = f 1+ ^L S j e -dT h /k_ 1 

< o 

and using l'Hospital's Rule 



lim (ke- dTh/k + dT h -k) = dT h + lim 

k — too \ / k—>ot 



= dTh + lim 

k— >oo 



dT± e -dT h /k 



= 

we deduce that sgn(ke~ dTh / k + dT h — k) > Therefore, 



sgn( ^ ) < (41) 



dT r 



□ 



We can deduce that we hit the smallest minimal value for the budget for the largest 
possible T r in this case that corresponds to when k = 1. This happens when the release 
frequency equals the partial harvest frequency. 



4.2 Discussion 

Figure [3] represents the analytical results obtained in the previous sections for a chosen set 
of parameters. The plot includes the two studied cases: either one of the partial harvest 
and the release period is an integer multiple of the other. 

Under this set of possible scenarios, increasing the frequency of release beyond the fre- 
quency of harvest requires that the total number of predators to invest in be higher than 
that when releases take place less or as often as partial harvests. In the latter case, to 
ensure pest eradication, the total budget of predators to invest in is fixed, independently of 
the release period. 

These results imply that it is clearly less costly to protect a greenhouse culture for lower 
frequencies of release. Of additional economic interest, in this case, the biological treatment 
is always combined with partial harvesting, so that there is little or no extra cost linked 
to the presence of workers on-site. However, we recall that [13] previously demonstrated 
that the higher the release frequency, the smaller the worst-case damages. Combining the 
results from both studies seems to indicate that the most profitable release strategy among 
the possibilities that have been considered is the one where releases are synchronized with 
the partial harvests. 
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0.4 



□ 
□ 
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Figure 3: Variation of the minimal number of predators required per budget year \i as a 
function of release to harvest period ratio. Parameters are given the values (in arbitrary 
units): a x — a y — 0.5, d = 1, and the rate of growth /'(0), functional response g'(0) and 
numerical response h'(0) with respect to time when the ecosystem is pest-free, i.e. x p (t) = 0, 
are all equal to 1. 



5 Conclusion 

The results obtained in this paper for the stability of the system are yet another confirmation 
that inundative control can be an effective means of suppressing low pest invasions in a 
greenhouse. This requires that a sufficient number of predators are introduced in the system 
as in, for instance, [3 [9]. 

Our study aimed to provide a control strategy in the protection of continuously grown 
crops that are partially harvested on a regular basis. We demonstrated that partial har- 
vesting had a non-negligible effect on biological control and needed to be taken into account 
when devising a control strategy in the case of such crops. 

We thus investigated the combined effects of releases and partial harvests in terms of 
the relative frequencies of their implementation. We considered the case where these two 
events occurred at periods such that one was the integer multiple of the other, and with the 
two events coinciding over the longer period. In particular, we found when releases were 
as frequent as or less than the partial harvests, the minimal budget did not depend on the 
period of release but instead on the harvest parameters, the growth function of the pest 
population, the mortality of the predators and the functional response. When releases were 
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more frequent than the partial harvests however, the minimal budget value increased with 
the increasing frequency of the releases, exceeding the constant value obtained for the less 
frequent case. Combined with the findings of [13] which pointed out that higher release 
frequencies led to the optimal control policy, we concluded that for the set of possibilities 
that was studied, the current best strategy is when release and harvest frequencies are equal. 

This approach has, however, its shortcomings. Since the integer multiple factor is key 
to calculating the minimal budget which would satisfy the stability conditions, it is not yet 
generalised to other scenarios where neither period is the integer multiple of the other. This 
would happen for instance at other rational non-integer ratios as irrational ones. It is highly 
likely that these intermediate ratios might induce other dynamics in the system. Whether 
they might stabilise it given even lower minimal budget values or favour chaos remains to 
be seen. Moreover it would be interesting to extend the results to the case where the two 
controls never coincide in spite of following a periodic pattern. This would be in the line of 
the work, for instance, of [H], where pesticide spraying - which is analogous to harvests - 
and releases are not synchronised. 

Nevertheless, we consider that our simplification already has its practical economical 
advantage. Indeed, coinciding periods imply little or no additional costs incurred in terms 
of labour: the task of predator release can be assigned to workers in charge of partial 
harvesting. Field-testing is now the next step required to validate the results of this paper. 
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Appendix 

Propositions Let Hypothesis]]] hold, then 

Proof: In order to compute the integral, we describe y p h(t) as j/ p ; l (jX'+)e~(*~ jTr ) in each 
time interval [iT r , (i + l)T r ], with y p h(iT+) given by © when y(nT^) = ?/*. This yields: 



Jo h y P h(t)dt =y^y ph (iT+) e~ d(t -^dt 

fe— 1 / t— 1 \ 7 1 

= E y* e ~ ldTr + ^ E e ^ dTr J ' e ^ 



i=0 \ 3=0 ' ■'" 



2/ l_ e -<wy T/UrZjtO 1- 



•(l- e -«T0 + ^(*-^S^) 



d 

-dT. 



K3* 



l-(l-"y)e- dT " 

(l-e-^)(l-a H )+a y fl-e"^) ) (l- 



(l- e -^)(l-(l-a y )e-^) 





r 1 




dT h \ 






|(l-(l-a H )e-^) 






" 1 y 




«y)e-* r *) 



/' 





^(l-(l-a y ) e - d ^ 




l(l-e- d ^) 


dk i 


f dT h \ 


l(l-(l-o y )e- 











" k j 






|(l-(l-a H )e-^) 






-3^ 


l(l-(l- 


% )e- dr h) 



a„e- £ ^)(l-e-^)+(fc(l-e-^))(l-(l-a !/ ) 



(l-e T i j(l-(l-Q: ! ,)e- dT fe) 



d 1 [ l-{l- ay )e-^H ) \ k ( x _ e -^ 



□ 
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Proposition 4 Let Hypothesis^ hold, then 



y pr (t)dt 



[iT r 1 - e- dTh 
~d~l- {l-u y )e- dT * 



Proof: In order to compute the integral, we describe y pr (t) as y pr (iT^)e ^ %Th ^ in each 
time interval [iTh, (i + l)Tk], with y pr (iT^) given by (fl2| when j/(mT r + ) = y* . This yields: 



A i -(i+l)T h 

e ~d(t-iT h ) dt 



Jo y P r(t)dt = ^ Vpr (iT+) I 

2=0 Jl '- - 

,_n JO 



pT r (l- e -^h) l-(l-q H ) fc e - fcdT h 

l-(l-a„) fc e- dT r d l-(l-a y )e- dTh 

fiT r (l-e- dT h) 
d(l-(l-a v )e- dTh ) 



□ 
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